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ABSTRACT 

We present spatially resolved stellar rotation velocity and velocity dispersion profiles form 
Keck/LRIS absorption-line spectra for 25 galaxies, mostly visually classified ellipticals, in three clus- 
ters at z ~ 0.5. We interpret the kinematical data and HST photometry using oblate axisymmetric 
two-integral /(-E, Lz) dynamical models based on the Jeans equations. This yields good fits, provided 
that the seeing and observational characteristics are carefully modeled. The fits yield for each galaxy 
the dynamical M/ L and a measure of the galaxy rotation rate. Paper II addresses the implied M/ L 
evolution. Here we study the rotation-rate evolution by comparison to a sample of local elliptical 
galaxies of similar present-day luminosity. The brightest galaxies in the sample all rotate too slowly 
to account for their flattening, as is also observed at z = 0. But the average rotation rate is higher at 
z « 0.5 than locally. This may be due to a higher fraction of misclassified SO galaxies (although this 
effect is insufficient to explain the observed strong evolution of the cluster SO fraction with rcdshift). 
Alternatively, dry mergers between early-type galaxies may have decreased the average rotation rate 
over time. It is unclear whether such mergers are numerous enough in clusters to explain the observed 
trend quantitatively. Disk-disk mergers may affect the comparison through the so-called "progenitor 
bias" , but this cannot explain the direction of the observed rotation-rate evolution. Additional samples 
are needed to constrain possible environmental dependencies and cosmic variance in galaxy rotation 
rates. Either way, studies of the internal stellar dynamics of distant galaxies provide a valuable new 
approach for exploring galaxy evolution. 

Subject headings: galaxies: clusters: individual (CL3C295, CL0016+16, CL1601-I-42) — galaxies: 
evolution — galaxies: formation — galaxies: kinematics and dynamics. 



L INTRODUCTION 

Elliptical galaxies form an important component of the 
galaxy population in the universe. In fact, they consti- 
tute the bulk of the population both at the high-mass 
end of the galaxy distribution and at the highest envi- 
ronmental densities (i.e., in clusters). An understanding 
of the formation and evolution of elliptical galaxies is 
therefore crucial for any successful theory of hierarchical 
structure formation in a cosmological context. 

The most direct way to constrain the evolution of ellip- 
tical galaxies is through observations of their photomet- 
ric and dynamical properties at intermediate and high 
redshifts. Several of the most successful techniques to 
use such data have relied on the statistical properties of 
globally averaged or characteristic galaxy quantities. For 
example, the fundamental plane is a correlation between 
the effective radius Toff, the average surface brightness 
/off inside the effective radius, and the velocity disper- 
sion a. The evolution of the fundamental plane can be 
used to place constraints on the formation redshifts of 
elliptical galaxies (as reviewed in van Dokkum & van der 
Marel 2006, hereafter vDvdMOG). Similarly, the color- 
magnitude relation is a correlation between the total 
magnitude and globally averaged color. Elliptical galax- 
ies form a "red sequence" that can also be used to con- 
strain the formation redshifts of elliptical galaxies (e.g., 
De Lucia et al. 2004; Mei et al. 2006). 



The success and popularity of studies based entirely 
on globally averaged or characteristic galaxy properties 
can be attributed at least in part to the relative ease 
with which such studies can be performed. In particular, 
no modeling or understanding of the internal structure 
of the sample galaxies is necessarily required. However, 
this is at the same time the main drawback of such stud- 
ies. To draw conclusions one is usually forced to make 
some implicit assumptions about the structure of the 
galaxy, or at least, that this structure does not evolve 
with time. Such assumptions remain untested until one 
actually progresses to more detailed studies of the re- 
solved galaxy structure. In fact, data that enable such 
studies are often available, but remain unused. For exam- 
ple, fundamental plane or color-magnitude relation stud- 
ies do not use available information on the axial ratios 
of individual galaxies, their resolved surface brightness 
or color profiles, or any spatially resolved information on 
the kinematical structure of the galaxy. 

One approach to learn more about the internal struc- 
ture of distant elliptical galaxies is to construct detailed 
dynamical models for high-quality spatially resolved pho- 
tometric and kinematic data. The required dynamical 
modeling tools exist, and high quality photometric data 
is available for many galaxies from the Hubble Space 
Telescope (HST) . The most challenging part is therefore 
to obtain the necessary kinematic data. However, this 
is quite possible with large aperture telescopes equipped 
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with multi-slit spectrographs (e.g., Kelson et al. 2000). 
Treu & Koopmans (2004) and Koopmans et al. (2006) 
obtained and modeled resolved kinematical profiles for 
some lensing galaxies. They constructed relatively sim- 
ple dynamical models (namely, parameterized spherical 
models) but this still allowed them to obtain several in- 
teresting new results on the dark halo component of these 
galaxies. There is no reason why similar methodologies 
or more detailed modeling approaches cannot be applied 
to other large samples of (non-lensing) early-type galax- 
ies. 

In vDvdM06 we presented spectroscopy with the Low 
Resolution Imager and Spectrograph (LRIS) on Keck of 
some two dozen galaxies in the intermediate-redshift (z « 
0.5) clusters CL3C295, CL0016+16 and CL1601+42. 
The sample galaxies were selected to be bright enough 
for spectroscopy, and visually classified from the HST 
images of Dressier et al. (1997) and Small et al. (1997) 
as early-type (and in most cases elliptical) galaxies. In 
vDvdM06 we only used the data for a study of the FP 
evolution of the three sample clusters. However, we ob- 
tained spectroscopic data that are deep enough to extract 
spatially resolved rotation curves and velocity dispersion 
profiles. Furthermore, in the multi-slit mask designs care 
was taken to tilt individual slits so that they were aligned 
preferentially with the major axes of the galaxies, inso- 
far possible given geometric restrictions. We present the 
kinematical profiles here, and we construct detailed dy- 
namical models to interpret them. 

Our models are axisymmetric with two-integral distri- 
bution functions of the form / — f{E,Lz). Here E is 
the energy and the angular momentum around the 
symmetry axis. The internal dynamics are calculated 
by solving the Jeans equations of hydrostatic equilib- 
rium. Models of this type have been constructed suc- 
cessfully for many nearby galaxies (e.g., Binney, Davies, 
& lUingworth 1990; van der Marel 1991; Magorrian et 
al. 1998; Corsini et al. 1999). Two-integral models are 

relatively simple in that they have vj^ = u^. So they 
can be thought of as the axisymmetric generalization of 
spherical isotropic models. Of course, more sophisticated 
modeling approaches do exist (e.g., those that built com- 
pletely general three-integral distribution functions using 
the technique of numerical orbit superposition; van der 
Marel et al. 1998; Gebhardt et al. 2003). However, the 
complexity of these models would be "overkill" as com- 
pared to the quality of the data that can be realistically 
obtained at intermediate redshifts. 

The models yield the rotation rate of each galaxy by 
fitting to the spatially resolved kinematical profiles. The 
results can therefore be used to study the evolution of 
the rotation rate of galaxies as a function of time. This 
evolution can yield new insights into the mechanisms by 
which the internal structure of elliptical galaxies gets 
modified through mergers. Mergers are believed to be 
characteristic of the hierarchical build up of all struc- 
ture in the Universe. However, the relative importance 
of different types of mergers remains poorly constrained. 
While the formation of elliptical galaxies through disk- 
disk mergers has received considerable attention histori- 
cally (e.g., Barnes 1988), much recent work has focused 
on the importance of "dry" mergers between early-type 
galaxies (van Dokkum 2005; Bell et al. 2006). Such merg- 



ers produce remnants with different rotational proper- 
ties. Disk-disk mergers tend to produce remnants with 
relatively rapid rotation (Naab & Burkert 2003; Naab, 
Jesseit & Burkert 2006b), whereas dry mergers tend to 
produce remnants with relatively little rotation (Naab, 
Khochfar & Burkert 2006a) . A study of the evolution of 
galaxy rotation properties can therefore contribute to an 
improved understanding of the actual mechanisms that 
have shaped todays elliptical galaxies. 

The inferred rotation rates of the sample galaxies 
also provide the potential to kinematically identify SO 
galaxies that have been visually misclassified as ellipti- 
cal galaxies. SO galaxies make up only 12% to 21% of 
the galaxy population in the three clusters studied here, 
compared to 50% to 60% for nearby clusters (Dressier et 
al. 1997). In the three clusters, visually-classified ellipti- 
cal galaxies outnumber SO galaxies by a ratio of ~ 3 : 1, 
whereas in nearby clusters SO galaxies outnumber ellipti- 
cal galaxies. Issues of misclassification have been raised 
as a potential contributing factor to this strong apparent 
evolution of the SO fraction (e.g., Fabricant et al. 2000). 
This may be relevant for CL3C295 and CL0016-fl6, 
given that their deficit of SO galaxies is compensated by 
an overabundance of elliptical galaxies (which constitute 
more than 50% of the population). In these two clus- 
ters the total early- type galaxy fraction (i.e., elliptical 
and SO galaxies combined) is only slightly lower than in 
nearby clusters. By contrast, in CL1601-I-42 the majority 
of galaxies are spirals, which outnumber the SO galaxies 
by a factor of 4-5. 

The layout of the present paper is as follows. Section [2] 
discusses the details of the modeling approach. Section[3] 
presents the results for the individual galaxies in the sam- 
ple. Section [4] discusses the inferred rotation properties 
of the intermediate-redshift sample galaxies, compares 
them to the rotation properties of early-type galaxies in 
the local Universe, and discusses the results. Section [5] 
presents a summary and discussion of the paper's main 
conclusions. 

The dynamical models also yield the mass-to-light ra- 
tios M/L for all the galaxies. These are used in van der 
Marel & van Dokkum (2006, hereafter Paper II) to study 
the M/L evolution of elliptical galaxies as a function of 
redshift. This complements and tests analyses based on 
Fundamental Plane evolution, such as those presented in 
vDvdM06. 

Throughout both papers we assume a cosmology with 
ilni = 0.27, flA = 0.73 (the values obtained by 
the Wilkinson Microwave Anisotropy Probe; Spergel et 
al. 2003) and Hq — 71 kms~^ Mpc~^ (the value obtained 
by the HST Cepheid Key Project; Section 7 of Freedman 
et al. 2001). 

2. DYNAMICAL MODELS 

For the two-integral dynamical modeling we have used 
the assumptions and software described previously by 
van der Marel et al. (1994b). The goal is to first find a 
light density distribution that is consistent with observed 
HST photometry, and then to find a dynamical solution 
of the Jeans equations that reproduces the observed Keck 
spectroscopy. Fitting the model to the kinematics yields 
both a normalized measure k of the galaxy's rotation 
rate and the dynamical M/L over the spectroscopically 
explored region. 
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The three-dimensional mass density is assumed to be 
oblate axisymmetric and of the form 



piR,z)^po{m/b)"[l + {m/b) 



2i<5 



z'Q'\ 

Here z is the symmetry axis of the galaxy and (i?, </>, z) 
are the usual cylindrical coordinates. All isodensity sur- 
faces are spheroids with axial ratio Q < I. The adopted 
mass density has no particular physical significance. It 
is merely a convenient parameterization in which the pa- 
rameters pq, b, a and d provide sufficient freedom to ade- 
quately fit the data available for the sample galaxies. The 
components of the force V$ (where $ is the gravitational 
potential) exerted by the matter with mass density given 
in equation ([1]) can be expressed as one-dimensional inte- 
grals (Binney & Tremaine 1987, eq. [2-88]) that are easily 
evaluated numerically. 

The Jeans equations for hydrostatic equilibrium in a 
two-integral model 
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For known p and V$ the equations can be solved numer- 
ically on a two-dimensional grid in the meridional plane. 
The amount of mean streaming in the model is uncon- 
strained by the condition of hydrostatic equilibrium. We 
therefore introduce a free parameter k that assigns an 
amount of rotational support according to 



.,111/2 



(3) 



For fc = the model is non-rotating, whereas for |fc| = 1 
the velocity dispersion tensor is isotropic and the model 
is a so-called "oblate isotropic rotator" . The maximum 
value of |fc| for a physical model is defined by the criterion 



that CT^ 
that 



is everywhere positive. This implies 
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For line-of-sight projection we adopt a Cartesian co- 
ordinate system (x,y,w) with w along the line of sight, 
{x, y) in the plane of the sky, and x along the major axis. 
The projected mass density is 



^{x,y) = / 

J — { 



p{x^ y, w) dw. 



(5) 



We assume that the mass-to-light ratio M/L (expressed 
throughout this paper in units of Mq/ Lq) is constant 
with radius in the galaxy. This is known to be a reason- 
able assumption when modeling kinematical data inside 
the effective radius of elliptical galaxies (van der Marel 
1991; Kronawitter et al. 2000). The projected light den- 
sity is then fi{x,y) — Yj[x,y) / {M / L). When viewed at 
inclination angle i, the projected contours are similar el- 
lipses with projected axial ratio Q', with 



cos^ i + sin 
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Two models of the same Q' and different inclination pre- 
dict exactly the same projected light density if they have 



the same values of 6, a, 5, andjoQ, where jo = pa/{M/L) 
is a luminosity density. 

The projected luminosity-weighted line-of-sight veloc- 
ity moments of order i are given by 
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In the present context we are interested only in the first 
two moments. For a two-integral system 



Vw = COS (j) sin i 
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The resulting line-of-sight integrals can be evaluated nu- 
merically once the Jeans equations have been solved. 

The observational setup is characterized by the slit 
width W, the pixel scale P in the spatial direction, and 
the seeing FWHM S. The model predictions {vl,){x,y) 
must therefore be convolved (in an intensity-weighted 
sense) both with a square tophat function of size W x P 
and with the seeing point-spread function (PSF). The 
latter was assumed to be Gaussian. For the convo- 
lutions we used a Monte-Garlo integration procedure. 
This ultimately yields for each spatial pixel (or set of 
binned pixels) along the slit the model predictions for 
the mean velocity V and for the second velocity moment 
= V'^ + a"^ , where a is the velocity dispersion. 
For a given inclination, the model parameters jo, b, a, 
S and Q can be fixed by fitting the projected light den- 
sity p{x, y) to the available galaxy imaging, which was 
calibrated to the rest-frame i3-band. To transform the 
best-fitting luminosity density of each model to units of 
Lq pc""^ we adopted a solar absolute B-band magnitude 
Mb,q = 5.48 (Binney & Merrifield 1998) and the relation 
between arcsec and pc dictated by the adopted redshift 
and cosmology. For the model thus obtained we solve the 
Jeans equations using the trial values M/L ~ 1 for the 
mass-to-light ratio and k = 1 for the rotational support. 
The results are projected along the line-of-sight and con- 
volved with the observational setup. Let the predictions 
thus obtained be V and ctrms ■ The predictions for other 
values of M/ L and k are then obtained by simple linear 

scaling: V ^ fc(M/L)i/2]/ and ^rms ^ {M/Ly/^anMS- 
Let the observed mean velocity and velocity dispersion 
be Vohs and (Jobs- The data- model comparison resid- 



sys 



k {M/Lf/^ V and 



[(M/L)a|Ms - (M/L)T/2]i/2. Here Fsys is 



uals are then 5V — Vohs 

(Jcr = CTobs - [(M/L)ct|, 

the systemic velocity of the galaxy. We minimize these 
residuals in a sense (i.e., squared, weighted with the 
formal uncertainties AVobs and Adobs, and summed over 
all data points) to obtain the best-fitting values of M/L, 
fc, and Vsys. The formal uncertainties on these quantities 
(i.e., those that reflect propagation of random uncertain- 
ties) follow from standard theory. If the formal best 
fit for |fc| exceeds the maximum value fcmax allowed for a 
physical model by equation ([4]), then we adopt the latter. 
In discussions of Vohs in the remainder of this paper it is 
to be understood that V^ys was always subtracted. The 
systemic velocities are of little interest by themselves (in 
the present context) and will not be discussed further. 

The only parameter that is not generally constrained 
by the data is the inclination, which determines the ax- 
ial ratio Q according to equation One possibility is 
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TABLE 1 

Intermediate-Redshift Galaxy Sample 
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Note. — Column (1) lists the galaxy ID and column (2) the morphological type, both from Small et al. (1997). Column (3) 
lists the redshift used for the modeling. For cluster members this is the redshift of the cluster from Dressier & Gunn (1992) 
and Dressier et al. (1999); for the field galaxy CL 1601-270 it is the observed redshift reported in vDvdM06. Column (4) 
lists the major axis position angle and column (5) the apparent ellipticity (e = 1 — Q'), both determined as described in 
Section [3.21 Column (6) lists the intrinsic axial ratio Q of our "standard" model, as defined by equation ([9|. Column (7) lists 
the corresponding inclination, calculated from equation Column (8) lists the seeing FWHM of the spectroscopic observations 
in arcsec, determined by modeling the observed galaxy intensity profiles along the slit as described in Section [3.31 For the four 
galaxies that were observed with two different slit masks we list only the best value of the seeing. Columns (9) and (10) list 
the velocity dispersion for an aperture with the size of the effective radius, and its formal uncertainty. Columns (11)-(15) list 
the quantities inferred from our dynamical models with the standard axial ratio and inclination. Columns (11) and (12) list 
the best-fitting mass-to-light ratio (in solar B-band units) and its formal uncertainty. Columns (13)-(15) list the best-fitting 

rotational support and its formal uncertainties If A+fc = then this indicates that the best fit is obtained when k is 

set to the maximum value that guarantees positivity of cr^, as given by equation (|4]). An asterisk in column (13) indicates that 
the galaxy was included in the rotational support analysis presented in Section |4l 



to assume that all galaxies are viewed edge-on, so that 
i = 90° and the intrinsic axial ratio is equal to the (ob- 
served) projected axial ratio: Q = Q' . Alternatively, one 
can assume that all galaxies have a small and fixed ax- 
ial ratio, Q — Qmin! we have taken Qmin — 0.4 which 
is approximately the smallest intrinsic axial ratio found 
for early- type galaxies (Tremblay & Merritt 1995). The 
edge-on and Q = Qmin models are useful in that they 
produce, respectively, the roundest and flattest models 
possible, so that they bracket the range of available pos- 
sibilities. However, since they both represent extremes, 
they are less useful for calculating the average properties 
of a large sample. For that purpose it is better to use 
models that are "correct" in an average statistical sense. 
To construct such models we used available knowledge 
on the projected axial ratio distribution f{Q') of nearby 
early-type galaxies. Tremblay & Merritt (1995) inverted 
this distribution to obtain the distribution of intrinsic ax- 
ial ratios f{Q)- They did so under various assumptions 
for the underlying geometry, and we use here their re- 
sults for an oblate axisymmetric geometry. Given /(Q), 



one can find^ the distribution f{Q\Q') that a galaxy of 
known projected axial ratio Q' has intrinsic axial ratio Q. 
In turn, this can be used to calculate the mean intrinsic 
axial ratio as function of projected axial ratio Q': 

{Q){Q') = I Qf{Q\Q')dQ. (9) 

In our models we set the intrinsic axial ratio Q of each 
galaxy to the value {Q){Q') appropriate for its observed 
projected axial ratio Q' . We calculate the corresponding 
inclination from equation ([6]). We will refer to the values 
thus obtained as the "standard" axial ratio and inclina- 
tion. Results discussed in the remainder of this paper 
are for these standard values. Our approach does not 
guarantee that we use the correct inclination for any in- 
dividual galaxy. However, it does imply that our models 

^ We calculated f{Q\Q') with a simple Monte-Carlo scheme in 
which for a very largo number of hypothetical galaxies we perform 
the following steps: (a) draw a random three-dimensional viewing 
direction; (b) draw an intrinsic axial ratio Q from the distribution 
f(Q)\ (c) calculate Q' from equation l(6)l; (d) add the results to a 
joint frequency distribution of the quantities Q and Q' . 
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Fig. 1. — Elliptical isophote fit results for the sample galaxies as a function of semi-major axis length. Shown are from left to right: 
isophotal surface brightness, transformed to rest-frame B-band and corrected for (I-I-2)* cosmological surface brightness dimming; isophotal 
ellipticity; and major-axis position angle. The curves in the surface brightness panels show the predictions of the best-fit mass density 
models of the form given by equation |(TJ. The horizontal lines in the ellipticity and position angle pane ls show the constant values used in 
the dynamical modeling (these are the weighted average values determined as described in Section l3.2l l. 



have the correct average intrinsic axial ratios when av- 
eraged over a large sample. Where relevant, we address 
the dependence of the results on inclination by referring 
to results obtained with the edge-on models or Q — Qmin 
models. 

3. RESULTS FOR INDIVIDUAL GALAXIES 
3.1. Sample 

For the modeling we used the sample of galaxies from 
vDvdM06. The galaxies reside in the clusters CL 3C295, 
CL 0016-1-1609, and CL 1601+4253, at redshifts z = 



0.456, 0.546, and 0.539, respectively (Dressier & Gunn 
1992; Dressier et al. 1999). The clusters were selected 
based on their visibility at the time of the Keck observa- 
tions, and because they are among the most SO deficient 
clusters in the MORPHS sample (Dressier et al. 1997). 
The MORPHS sample itself was not selected according 
to strict criteria. The galaxy selection was largely con- 
strained by the geometry of the Keck/LRIS masks, and 
by the fact that sample galaxies should be bright enough 
for spectroscopy. Priority was given to galaxies classified 
from HST images as E or E/SO by Small et al. (1999). 
The latest-type galaxy included in the sample was an 
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Fig. 1. — (continued) 



SO/Sb galaxy. This galaxy 3C295-568 was included for 
the specific purpose to see if rotation could be measured 
reliably. For the present paper we excluded two galaxies 
from the sample of vDvdMOG: 3C295-47 and CL1601- 
474. The former is a foreground field galaxy at z = 0.131. 
The latter has a disturbed morphology suggesting an on- 
going merger, so that equilibrium models are likely to 
be of limited value. One other galaxy in the sample of 
vDvdM06 was found to be a field galaxy, namely CL 
1601-270. This galaxy was retained because its redshift 
z = 0.510 is close to that of the other galaxies. 

The galaxies in the final sample are listed in Table [H 
including redshifts, visual classifications, and aes- The 
effective velocity dispersions are included for use in Pa- 



per II. They were calculated from the observed veloc- 
ity dispersions measured from spectra that summed the 
central five pixels along each slit (presented in vDvdM). 
These measurements were transformed to an estimate of 
the dispersion inside an aperture of size r^s using the 
formulae given in Jorgensen et al. (1995b). The values 
of Toff and /off for the sample galaxies can be found in 
vDvdMOe. 

3.2. Photometry 

To start the analysis we used the HST/WFPC2 images 
previously shown and discussed in vDvdMOG. For the 
galaxies in the clusters CL3C295 and CL160H-42 there 
are only images in the F702W filter. For the galaxies in 
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Fig. 1. — (continued) 



-1 -.5 
log(r[arcscc]) 



.6 — 
A — 
.2 




-.6 -A -,3 
log(r[arcscc]) 



100 
4 50 



-.6 -A -.3 
log(r[arcscc]) 



the cluster CL0016+16 there are images in both F555W 
and F814W, but we only used the latter. The images 
were corrected for the blurring of the HST PSF using 
Lucy-Richardson deconvolution (Richardson 1972; Lucy 
1974). The PSF used for the deconvolution was calcu- 
lated with the TinyTim software package (Krist & Hook 
2001). The deconvolved galaxy images were analyzed 
through isophotal ellipse fitting with the task ELLIPSE 
in the IRAF software package (Jedrzejewski 1987). This 
yields for each galaxy the major axis profile of surface 
brightness, ellipticity, and position angle. The results 
are shown in Figure [TJ The profiles can generally be 
followed out to 1-5 arcsec from the galaxy center. The 
ellipticity and position angle are not shown in the figure 



inside 0.1 arcsec, where they are generally unreliable due 
to the residual effects of the HST PSF. They are also not 
shown outside of 1.0 arcsec, where they tend to become 
unreliable due to the low galaxy brightness relative to 
the sky brightness. 

The models described in Section [5] are axisymmetric 
with constant axial ratio, and in projection have con- 
stant ellipticity e = 1 — Q' and major axis position angle 
PAmaj- To optimally choose these parameters for our 
models we calculated for each galaxy the average ellip- 
ticity and major axis position angle between 0.1 and 1 
arcsec, weighted both with the flux and the annular area 
that each isophote represents. The results are listed in 
Table[T]and are also shown as horizontal lines in Figure[TJ 
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Fig. 2. — Kinematics of the sample galaxies. Each set of two adjacent panels show the rotation velocity V and the velocity dispersion 
tr along a given slit. The galaxy name is indicated in the label above the rotation velocity panel. The slit position angle is indicated in 
parentheses and is expressed as the sum of the major axis position angle given in Table [T] and an offset angle. The abscissa is given in 
arcsec. Data points are indicated with error bars. The heavy solid curves show the predictions of the best-fitting models with the standard 
inclination (see Table [TJ, taking into account the observational setup and spatial resolution. For comparison, the dashed curves show the 
model predictions as they would appear at infinite spatial resolution. 



Clearly, the assumption of constant ellipticity and posi- 
tion angle with radius is better for some galaxies (e.g., 
CL1601-2040) than for others (e.g., CL0016-2050). How- 
ever, even for galaxies with significant ellipticity varia- 
tions or isophotal twists our approach still has the ad- 
vantage that it properly takes into account some average 
measure of galaxy flattening. 

The observed surface brightnesses in counts per sec- 
ond per pixel were transformed to magnitudes per square 
arcsecond using the known WFPC2 pixel sizes and filter 
zero-points from the WFPC2 Data Handbook (Baggett 



et al. 2002). The results were corrected for foreground 
Galactic extinction using the Schlegel et al. (1998) maps. 
The magnitudes in the filters used for the observations 
were then transformed to galaxy rest-frame B-band mag- 
nitudes using the relations given in vDvdM06. And fi- 
nally, a (1-l-z)'' correction for cosmological surface bright- 
ness dimming was applied. 

For each galaxy we determined the parameters so that 
the projected luminosity density fi(x, y) of the model 
defined by equation ([T]) best fits the observed surface 
brightness profile. The fit was done using a Levenberg- 
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Fig. 2. — (continued) 



Marquardt iteration scheme (Press et al. 1992). Data 
points within 2 arcsec from the galaxy center were in- 
cluded in the fit. The integrated brightness of the model 
in the central pixel was also included. The predictions 
of the best model fits are overplotted on the data in Fig- 
ure [TJ The predictions generally fit the observed surface 
brightness profiles reasonably well; the median residual 
over the entire sample is 0.09 mag/arcsec^. The worst 
fit is obtained for CL 3C295-834, which clearly has more 
structure in its profile than can be accounted for by our 
model. The adopted inclination i and axial ratio Q for 
our standard model of each galaxy are listed in Table [TJ 
For a given galaxy, the luminosity density jo is inversely 
proportional to Q because the total luminosity of the 



model, which is fixed by the observations, is proportional 
to joQ. However, the inferred parameters &, a, and S (see 
eq. [T]) and the fits shown in Figure [1] are all independent 
of the assumed i and Q. 

The galaxy CL 3C295-2014 is the well-known AGN 
3C295. There is the possibility that this galaxy contains 
a central non-thermal point source that could bias the 
analysis. However, the results for both the light profile 
of this galaxy and its M/L (see Figure[T]and Paper II) do 
not provide any evidence for deviations from the trends 
defined by the other galaxies in the sample. We therefore 
retained CL 3C295-2014 in our sample and did not treat 
it in any special way. 
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Fig. 2. — (continued) 



3.3. Kinematics 

For each galaxy we solved the Jeans equations to ob- 
tain predictions for the projected kinematics on the plane 
of the sky. To model the observational setup we used 
the known slit width = 1.1 arcsec and pixel scale 
P = 0.215 arcsec/pixel. Binning along the slit that was 
applied to the data (see vDvdM06) was also taken into 
account. We assume that during the observations all slits 
were perfectly centered on their corresponding galaxy. 
This was probably true to within ~ 0.2" accuracy. To 
determine the seeing for the observations with a given slit 
mask we ran models with different values of the FWHM 
S. We compared the observed and predicted intensity 
profiles along the slit for the galaxies observed with that 
mask, and adopted the value that gave the best agree- 
ment in a sense. The inferred seeing FWHM values 
range from 0.71" to 0.91" (consistent with estimates ob- 
tained from inspection of stars observed during the same 
run) and are listed in Table [H The position angle of 
each slit with respect to the galaxy major axis was deter- 
mined by subtracting the value PA„iaj in Table[T]from the 
known slit position angle on the sky. The values of the 



mass-to-light ratio M/L and rotational support k were 
determined to optimize the fit to the observed kinematics 
along the slit. For four galaxies (namely: CL0016-461, 
CL1601-524, CL1601-753, CL1601-814) data were avail- 
able for more than one slit position angle (from observa- 
tions with different slit masks). In those cases a single 
combination of M/L and k was fit simultaneously to the 
data for both slits. 

Figure [2] shows the resolved kinematical profiles for all 
galaxies (obtained as described in Section 2.6 of vD- 
vdM06), together with the best fits of the standard- 
inclination models. The fits are generally good, as con- 
firmed by the values. When summed over the en- 
tire sample, the mean velocity data points are fit with 
= 145.8 for 113 degrees of freedom. The velocity dis- 
persion data points are fit with = 81.1 for 102 degrees 
of freedom. Models constructed with other assumptions 
for the inclination are both visually and statistically in- 
distinguishable from those constructed for the standard 
inclinations. Comparison of edge-on models and models 
with Q = Qmin yields a difference in total summed 
over the entire sample, of only Ax^ = 1.0. This indi- 
cates that our kinematical data provide insufficient in- 
formation to constrain the galaxy inclinations. Indeed, 
that would require high quality two-dimensional velocity 
fields (Cappellari et al. 2006). 

For comparison. Figure [2] also shows the model predic- 
tions before convolution with the slit width, pixel size, 
and PSF. It is clear that this convolution has a major im- 
pact on the predictions, due to the fact that the spatial 
resolution of the observations is not much smaller than 
the galaxies themselves. In particular: (a) the observed 
rotation curves are always much shallower than the rota- 
tion curves that would be observed with infinite spatial 
resolution; and (b) the velocity dispersion profiles fall off 
less steeply with radius than they would when observed 
with infinite spatial resolution. Both effects are easily un- 
derstood as the consequence of scattering of light from 
small radii (where the galaxy is brightest) to larger radii. 
Obviously, the resolved kinematics of distant galaxies can 
only be meaningfully interpreted if detailed dynamical 
models are combined with a proper accounting of the 
spatial resolution of the observations, as we have done 
here. 

The kinematical data typically extend to between 0.2" 
and 1.2" from the galaxy center with a typical value of 
^ 0.7". At z = 0.5 this corresponds to 4.3 kpc. So 
our data extend about as far out as the data available 
for many nearby elliptical galaxies (e.g., van der Marel 
1991; Emsellem et al. 2004). On the other hand, the sam- 
pling and resolution of the kinematical profiles in physi- 
cal units are of course much worse. Nonetheless, for the 
mean velocities V it is clear that our spatially resolved 
measurements do provide important new insight on the 
internal structure of the galaxies as compared to a single 
integrated spectrum. The results for the SO/Sb galaxy 
3C295-568 show that rotation can be reliable measured. 
While this is the most rapidly rotating galaxy in the sam- 
ple (as expected, given its visual classification), most of 
the other galaxies also have statistically significant gra- 
dients in V along the slit. With the help of our models 
this provides constraints on the intrinsic rotation veloc- 
ity. By contrast, for the velocity dispersions a it is not 
clear that we are able to learn much more about the in- 
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Fig. 3. — (a) (vja)* versus B-band absolute magnitude Mb for nearby elliptical galaxies. The data are from the compilations of Davies 
et al. (1983) and Bender, Burstein &; Faber (1992). The values from these papers were transformed to Hq = 71 kms~^ Mpc~^, for 
consistency with the value assumed here, (b) k versus Mb for the sample of intermediate redshift cluster galaxies. The values of Mg were 
corrected to z = using the luminosity evolution calculated by vDvdM06, for comparison with the local sample. Symbol types indicate the 
visual galaxy classification from Smail et al. (1997): E (solid dot), E/SO (hexagon), SO/E (square), SO (triangle), or SO/Sa (open circle). 
Galaxies with apparent ellipticity e < 0.10 were excluded from both panels. For the intermediate redshift cluster sample, galaxies were 
included only if there is data along a slit within 45° from the major axis. In both panels, galaxies with Mg < —20.4 span a range of 
(v/a)* or k values. However, most of these high-luminosity galaxies do not rotate fast enough to account for their flattening (i.e., they 
have (v/a)* < 1 or fe < 1). 



ternal structure than we would from a single integrated 
measurement. The model predictions in Figure [5] show 
that any gradients in a over the central arcsec are ex- 
pected to be small. To discriminate between different 
models one would need to be able to reliably measure 
these small gradients. In most galaxies our formal un- 
certainties on a are not small enough to allow this. Of 
course, we do measure an accurate mass-to-light ratio, 
but that would have been possible also with only a sin- 
gle integrated a measurement. 

Table [1] lists the best- fitting values of k and M/L 
(in rest-frame i?-band solar units) for the standard- 
inclination models, with their formal uncertainties. The 
sign of k is of little interest and indicates merely in which 
of two possible directions the galaxy is rotating around 
its symmetry axis. Adding 180° to the value of PAmaj 
yields a reversal of the sign of k. For convenience we have 
defined the values of PAmaj in Table[T]to be such that the 
best-fitting value of k is always positive. The M/L val- 
ues of the galaxies are not addressed further here, since 
they are discussed and interpreted in detail in Paper II. 
Note that M/L depends on the assumed distance, and 
therefore on the assumed cosmological parameters flm, 
17 A and Hq (as listed in Section [T]) . The uncertainties 
in M/L introduced by distance uncertainties are not in- 
cluded in Table [T] but they are quantified in Paper II. 

4. GALAXY ROTATION PROPERTIES 

4.1. Intermediate- Redshift Cluster Sample 

The model quantity fc is a normalized measure of the 
rotation rate of a galaxy, with an oblate isotropic rotator 
model corresponding to A; = 1. Inspection of Table [1] 
shows that the random uncertainties in k tend to be large 
when it is estimated from data along a slit that is far from 
the major axis, or for any galaxy that is nearly circular 
in projection on the sky. In these cases the value of k 
is also very sensitive to small systematic uncertainties in 



the measured major axis position angle or the ellipticity 
of the galaxy. For this reason we consider in the following 
only the 15 galaxies in the sample of intermediate redshift 
cluster galaxies for which a slit was placed within 45° 
of the major axis, and which have projected ellipticity 
e > 0.10. These galaxies are marked with an asterisk in 
column (13) of Table [T] 

Figure [Sb shows the inferred values of fc as a function 
of the B-band absolute magnitude Mb ■ The latter was 
obtained from the total B-band luminosity, which was 
estimated as Lb = 27rr^jj/off . For comparison to local 
samples, the luminosity was corrected to ^ = using 
the luminosity evolution calculated by vDvdM06. The 
galaxies span a range of fc values between to about 1.3. 
The galaxies that were not visually classified as E galax- 
ies are indicated with special symbols in the figure. The 
galaxies CL 0016-650 (E/SO), CL 1601-270 (SO/E) and 
CL 3C295-868 (SO) all have only modest rotation rates 
of fc = 0.3-0.4, albeit with sizable uncertainties of similar 
magnitude. Either way, the rotation rates of these galax- 
ies do not make them stand out among the ellipticals. By 
contrast, the galaxy CL 3C 295-568 (SO/Sb) has well- 
measured (Afc < 0.01) rapid rotation (fc w 1.08). This 
galaxy also distinguishes itself among the sample galaxies 
by having the highest apparent flattening (e = 0.58). 

4.2. Comparison to Local Elliptical Galaxies 

The quantity that is most often measured for local 
galaxies to quantify rotational support is (v/a)*. This 
is defined as the ratio of the the observed (v/a) and the 
value of {v/a) expected for an oblate isotropic rotator 
model with the observed apparent fiattening. This is al- 
most identical to the model parameter fc used here. The 
main difference is that {v/a)* is defined in terms of pro- 
jected global quantities (maximum rotation velocity and 
average velocity dispersion), whereas fc is defined locally 
by equation In principle this difference in defini- 

tion could lead to small systematic differences between 
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k and (v/a)*, even when measured on the same galaxy. 
However, these differences are unlikely to exceed '^0.1. 
This is negligible compared to the typical uncertainties 
in our k measurements for the intermediate redshift clus- 
ter galaxies (see FigurelSb). In the following we therefore 
compare the quantity k directly to the quantity (v/a)* 
available for local galaxies. 

As local comparison sample we combined the litera- 
ture compilations of {v/a)* values for elliptical galaxies 
from Davies et al. (1983) and Bender, Burstein & Faber 
(1992). The resulting sample contains 80 galaxies. For 
consistency with our treatment of the sample of inter- 
mediate redshift cluster galaxies we excluded from the 
following discussion the 12 galaxies with e < 0.10. 

Figure [3^ shows (v/a)* for the local ellipticals as a 
function of Mb- The local galaxies show a well-known 
trend between (v/a)* and Mb- Low-luminosity ellipti- 
cals {Mb > —20.0) have values of {v/a)* that cluster 
near unity. This indicates that these galaxies are typi- 
cally supported by rotation. High-luminosity ellipticals 
{Mb ^ —20.0) have values of {v/a)* that range from zero 
to just above unity, but with most of the galaxies clus- 
tering near zero. This indicates that these galaxies are 
typically supported by velocity dispersion anisotropy (for 
a more detailed discussion of what {v/a)* measurements 
imply for the velocity dispersion anisotropy of elliptical 
galaxies see Cappellari et al. 2005). The three galaxies 
in Figure [3^ with k w 1.5 and Mb between —20.0 and 
-21.0 are NGC 584, 2974, and 3640. These are proba- 
bly not normal ellipticals. NGC 584 is actually classified 
in the RSA catalog (Sandage & Tammann 1981) as an 
SO galaxy, although it is an E4 in the RC3. NGC 2974 
is also an E4 galaxy, but it has all the hallmarks of be- 
ing a misclassified SO galaxy (Cinzano & van der Marel 
1994). And NGC 3640 has considerable fine-structure 
that suggests it is a recent merger remnant (Prugniel et 
al. 1988). 

The galaxies in the sample of intermediate redshift 
cluster galaxies are all brighter than Mb ~ —20.0. The 
fact that they span a wide range of k values is therefore 
consistent with the behavior seen in the local sample. 
Moreover, the four brightest galaxies {Mb < —22.0) all 
appear to have insufficient rotation to account for their 
flattening, consistent with the situation for the elliptical 
galaxies in the nearby universe. However, the details of 
the distribution of rotational support are somewhat dif- 
ferent between the local and distant samples. Figure |4] 
compares a histogram of the k values for the intermedi- 
ate redshift cluster galaxies to a histogram of the {v/a)* 
values of local galaxies. For the intermediate redshift 
cluster galaxies we included only the 11 galaxies clas- 
sified as ellipticals, for consistency with the local sam- 
ple. For the local galaxies we included only galaxies with 
Mb < —20.4. With this selection, both histograms probe 
similar absolute magnitude ranges (for the local sample 
Mb has an average of —21.2 with a dispersion of 0.7 and 
for the distant sample it has an average of —21.5 with a 
dispersion of 0.9). 

The histograms in Figure|3]show that in a relative sense 
there is a higher number of rapidly rotating galaxies in 
the distant sample than in the local sample. The average 
(fc) for the distant galaxies is 0.68 ±0.13, whereas the av- 
erage {{v/a)*) for the local galaxies is 0.47±0.07. In the 
distant sample 36% of the galaxies have k > 0.9, whereas 




k or {v/a)' 

Fig. 4. — Comparison of a histogram of the k values for the in- 
termediate redshift cluster galaxies (heavy line) to a histogram of 
the (v/a)* values of local galaxies (dotted line). For the intermedi- 
ate redshift cluster galaxies we included only the 11 galaxies with 
reliable k determinations that were classified as ellipticals, for con- 
sistency with the local sample. For the local galaxies we included 
only galaxies with Mg < —20.4, so that both histograms probe 
similar absolute magnitude ranges. The histograms show that in a 
relative sense there is a higher number of rapidly rotating galaxies 
in the distant sample than in the local sample. 

in the local sample only 19% have {v/a)* > 0.9. This es- 
timate of the percentage difference at high rotation rates 
is actually conservative for two reasons. First, the per- 
centage of the local sample with {v/a)* > 0.9 would drop 
by 5% if we excluded the two potential SO galaxies NGC 
584 and 2974. And second, values of {v/a)* are defined 
as the rotational support for an edge-on model. In con- 
trast, our standard two-integral models are somewhat 
inclined. If instead we had used the values of k inferred 
from edge-on models for the intermediate redshift clus- 
ters, then all k values would have gone up by ~ 7%. The 
discrepancy between the histograms in Figure [4] would 
then have been even larger. In view of the statistics and 
these additional considerations, we conclude that the dif- 
ference between the rotation properties of the galaxies in 
our z « 0.5 cluster sample and those of our local compar- 
ison sample are significant. In the following subsections 
we discuss possible explanations for this difference. 

4.3. SO Mis classification 

One plausible explanation for the relatively high frac- 
tion of rapidly rotating galaxies at z « 0.5 is that the 
fraction of misclassified SO galaxies among visually clas- 
sified elliptical galaxies is larger in intermediate-redshift 
clusters than it is locally. This has in fact been suggested 
(e.g., Fabricant et al. 2000) as a possible explanation for 
the low fraction of SO galaxies found in intermediate- 
redshift clusters (Dressier et al. 1997). The bulges of 
SO and spiral galaxies typically have rotation rates con- 
sistent with oblate isotropic rotator models (Davies et 
al. 1983). Any contamination from disk light will further 
increase the observed rotation. Therefore, one would nor- 
mally expect to find SO galaxies at values of k or {v/a)* 
near or above unity. 

Based on the statistics for the local comparison sam- 
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pie one would have expected 0.19 x 11 = 2.1 galaxies in 
the distant galaxy histogram in Figure [4] with k > 0.9. 
Instead, 4 out of 11 are observed to have such high ro- 
tation rates. We excluded 4 galaxies with e < 0.10 from 
the histogram. These are unlikely to be SO galaxies, 
since the solid angle under which galaxies are viewed 
pole-on is small. Hence, we estimate that 13% (namely, 
(4 — 2.1)/(ll-|-4)) of the elliptical galaxies in these inter- 
mediate redshift clusters might be misclassified SO galax- 
ies. The average of the morphological percentages re- 
ported by Dressier et al. (1997) for the three clusters 
studied here are E:SO:Spirals = 47:16:37. Upon correc- 
tion for the suggested misclassification of SOs as ellipti- 
cals we obtain instead 41:22:37. Although correcting for 
misclassifications increases the SO fraction of these clus- 
ters, it remains much smaller than the ~ 60% found in 
local clusters (Dressier 1980). So there appears to be 
a strong evolution in the real fraction of SO galaxies in 
clusters. Our data do not suggest that this can be at- 
tributed to misclassification, although it should be kept 
in mind that the SO fraction in intermediate redshift clus- 
ters might be higher than what we just estimated. Not 
all SO galaxies would necessarily be recognized from their 
kinematics alone. 

Our results with respect to potential SO misclassifi- 
cation probably depend on the sample selection. The 
three clusters in our sample were selected based on their 
low fractions of visually-classified SO galaxies (see Sec- 
tion 13. ip . On the one hand, this makes these clusters 
good test cases: there are many visually-classified ellip- 
tical galaxies that could potentially be misclassified SO 
galaxies. On the other hand, if these clusters are truly 
poor in SO galaxies, then they are not good test cases: 
in that case there simply aren't many SO galaxies that 
one could potentially misclassify. Therefore, the results 
presented here may not generalize to other clusters at 
similar redshift that have higher SO fractions based on 
visual classification. 

4.4. Redshift Evolution of Galaxy Rotational Support 
through Mergers 

An alternative interpretation for the difference between 
the histograms in Figure |4] is to assume that there is a 
true redshift evolution in the internal dynamics of el- 
liptical galaxies. Since elliptical galaxies are collision- 
less systems, mergers are required to significantly change 
their rotation rate. There are different types of mergers 
and these can affect the evolution of the observed rota- 
tion rate in different ways. We will discuss the relevant 
types in turn. Independent of the type of merger, we are 
only concerned here with mergers that happen between 
z w 0.5 and z = 0, i.e., in the last 5 Gyr. 

Mergers between spiral galaxies form remnants that 
resemble elliptical galaxies (e.g., Barnes 1988). The pro- 
genitors of such mergers are not included in samples of 
elliptical galaxies at intermediate redshift. However, the 
remnants of such mergers are included in local samples 
of elliptical galaxies. Such mergers therefore contribute 
to the so-called "progenitor bias" (e.g., van Dokkum & 
Franx 2001). Dissipationless and dissipational simula- 
tions show that the remnants of disk-disk mergers have 
properties that are similar to local low-luminosity ellip- 
ticals, which tend to be disky and flattened by rotation 
(Naab & Burkert 2003; Naab et al. 2006b). In con- 



trast, disk-disk mergers do not seem a plausible forma- 
tion mechanism for the most luminous ellipticals, since 
only equal-mass mergers with special initial orientations 
can produce purely boxy anisotropic merger remnants. 
So if it were possible to remove from the local sample in 
Figure |4] the ellipticals that formed from disk-disk merg- 
ers since z « 0.5, then this would probably decrease the 
relative proportion of rapidly-rotating ellipticals. There- 
fore, progenitor bias due to disk-disk mergers cannot ex- 
plain the evolution seen in Figure [H In fact, the progen- 
itor bias is such that the observed evolution is probably 
only a lower limit to the true evolution. 

The situation is different for mergers between elliptical 
galaxies. The progenitors of such mergers are included 
in samples of elliptical galaxies at intermediate redshift. 
Therefore, such mergers introduce no progenitor bias 
(provided that luminosity-dependent sample selection ef- 
fects are properly understood). Simulations show that 
the remnants of mergers between early-type galaxies have 
properties that are similar to luminous ellipticals in the 
local universe (Naab et al. 2006a). In particular, the ro- 
tation rate of the remnants is generally low, even when 
the progenitors were flattened by rotation. Because el- 
liptical galaxies have little gas, there is no formation of a 
gas disk in the merger remnant that can form stars and 
boost the observed rotation rate (in contrast to the situa- 
tion for disk-disk mergers; Naab et al. 2006b). Therefore, 
disspationless ( "dry" ) mergers between elliptical galaxies 
can produce evolution in the elliptical galaxy population 
towards lower rotation rates. This would be consistent 
with the evolution seen in Figure U) 

The importance of mergers for the evolution seen in 
Figure [4] depends on the number of mergers of the rele- 
vant types in the last 5 Gyr. There is observational evi- 
dence that merging is indeed happening at an important 
rate in this time interval. Most studies originally lumped 
various types of interactions and mergers together (e.g., 
Le Fevre et al. 2000; Patton et al. 2002; Conselice et 
al. 2003; Lin et al. 2004; Bundy et al. 2004). How- 
ever, some recent studies have aimed to separate merg- 
ers of different types to specifically demonstrate the im- 
portance of dry merging (e.g., van Dokkum 2005). Bell 
et al. (2006) used data for galaxies in the GEMS sur- 
vey to estimate that present-day spheroidal galaxies with 
My < —20.5 have on average undergone between 0.5 and 
2 dry mergers since z « 0.7. Semi- analytical models of 
hierarchical structure formation predict that in the mass- 
range of interest to the present study, disk-disk mergers 
are about as common as dry mergers (Naab et al. 2006a) . 
Both types of mergers are therefore likely to have been 
numerous enough in the last 5 Gyr to have affected a 
significant fraction of elliptical galaxies in the luminosity 
range covered by our sample. 

Most observational work of the merger rate as a func- 
tion of redshift has studied field galaxies. While these 
studies have shown that mergers are important, this does 
not necessarily imply that our results have been affected 
by them. This is because the frequency of merging de- 
pends on environment. Our local comparison sample of 
elliptical galaxy rotation rates (see Section 14. 2|) is com- 
posed mostly of galaxies in the field and in low-density 
cluster environments such as Virgo. This range of envi- 
ronments is probably not too dissimilar from that probed 
by field galaxy surveys out to significant redshifts. This 
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assumption implies that disk-disk mergers in the last 5 
Gyr have probably been numerous enough among our 
local comparison sample to have affected the results in 
Figure [4] through the progenitor bias previously eluded 
to. 

Whether dry mergers have been numerous enough to 
decrease the rotation rates of the elliptical galaxies in our 
intermediate-redshift sample is a different issue. This 
is because those galaxies reside in dense clusters. The 
traditional belief has been that galaxy interactions in 
clusters are different from those in the field due to the 
large relative velocities between galaxies. Galaxy trans- 
formations can then occur due to the cumulative effect 
of frequent high speed galaxy encounters ("galaxy ha- 
rassment"; Moore et al. 1996). However, calculations of 
merger rates (based mostly on somewhat idealized sim- 
ulations without a cosmological context; e.g., Makino & 
Hut 1997) show that true mergers may be rare (except for 
accretion events onto the brightest cluster galaxy). How- 
ever, recent observations have challenged this paradigm. 
In particular, van Dokkum et al. (1999) and Tran et 
al. (2005) found direct observational evidence for ongoing 
mergers in the cluster MS 1054-03 at z = 0.83. The true 
abundance of dry mergers in clusters therefore remains 
unclear. However, it does appear possible that dry merg- 
ers in clusters may have caused a decline in the rotation 
rate of elliptical galaxies over the last 5 Gyr. 

4.5. Environmental Dependence and Cosmic Variance 

A final possible explanation for the difference between 
the rotation rates of our sample galaxies at z « 0.5 and 
those of the local comparison sample of Section HT^ is that 
the average rotation rate of elliptical galaxies may not 
be the same at different locations in the Universe. There 
could be a dependence on environment or there could 
be significant statistical fluctuations between widely sep- 
arated locations (cosmic variance). Our intermediate- 
redshift sample consists of galaxies in dense clusters. In 
contrast, our local comparison sample contains an in- 
homogeneous mixture of galaxies in cluster, group and 
field environments. A dependence on environment would 
therefore be particularly relevant for the comparison. 
However, observations in the local universe have pro- 
vided no evidence that such a dependence exists. The 
(v/a)* distributions of the rare elliptical galaxies in very 
low-density environments are quite similar to those for 
local samples biased towards high-density environments 
(Hau & Forbes 2006). Current cosmological models do 
not have the resolution or detailed physics that would 
be needed to make confident predictions for the size of 
either an environmental dependence or cosmic variance 
in the elliptical galaxy rotation rate. A good way to ad- 
dress this will therefore be to study larger samples of 
intermediate-redshift galaxies in a wider range of loca- 
tions and environments. 

5. SUMMARY AND DISCUSSION 

We have presented detailed dynamical models for 
the spatially resolved kinematics of 25 galaxies in 
the intermediate-redshift {z « 0.5) clusters CL3C295, 
CL0016-I-16 and CL160H-42. The sample galaxies were 
selected to be bright enough for spectroscopy, and vi- 
sually classified from HST images as early type (and in 
most cases elliptical) galaxies. The spectroscopic data 



were obtained with LRIS on Keck as described in vD- 
vdM06, with a 1.1" wide slit in seeing of 0.71" to 0.91" 
FWHM. The data sample a region that extends typically 
to ~ 0.7" from the galaxy center 4.3 kpc) sampled at 
0.215" per pixel. 

To model the data we adopted an oblate axisymmet- 
ric luminosity density of fixed axial ratio for each galaxy 
that fits the surface brightness profile inferred from HST 
imaging. The inclination of the models was chosen to 
provide the most likely intrinsic axial ratios in a statis- 
tical sense. Edge-on models and highly-flattened mod- 
els were also explored, but the results do not depend 
strongly on the assumed inclination. For each galaxy 
we solved the Jeans equations of hydrostatic equilibrium 
under the assumption that the distribution function has 
the two-integral form / = f{E, L^)- The resulting model 
predictions were projected along the line of sight and 
convolved with the observational seeing, slit width, and 
binning along the slit. Accounting for the observational 
characteristics is crucial for a proper interpretation of re- 
solved kinematical data obtained at these redshifts, and 
in particular the rotation curves. When this is done, the 
model predictions overall provide good fits to the kine- 
matical data. 

Comparison of models and observations yields for each 
galaxy a quantity k that is a normalized measure of the 
rotation rate of the galaxy. This quantity is similar to 
the quantity (v/a)* that is often used to measure the ro- 
tation rate of nearby galaxies (where k = {v/a)* = 1 for 
an oblate isotropic rotator model). To quantify the rota- 
tional properties of the sample galaxies we studied how 
k behaves as a function of galaxy magnitude Mb- We 
compared this to the known correlations between {v/a)* 
and Mb for elliptical galaxies in the local Universe. We 
find that the intermediate-redshift ellipticals span a large 
range in k from about to 1.3. This is not dissimilar from 
what is seen locally for elliptical galaxies of similar lumi- 
nosity. The four brightest galaxies in the sample, which 
will have Mb < — 22 even when evolved forward in time 
to z = 0, all rotate too slowly to account for their flatten- 
ing. This is consistent with what is seen in the local Uni- 
verse for galaxies of the same luminosity. Nonetheless, 
in a relative sense there is a higher number of rapidly ro- 
tating galaxies in the intermediate-redshift sample than 
in the local Universe (when comparing galaxies of similar 
present-day luminosity) . 

The general similarity of the rotation rates of the sam- 
ple galaxies to those of local ellipticals suggests that 
the intermediate-redshift galaxies have in large major- 
ity been correctly classified visually as early-type galax- 
ies. Nonetheless, a possible explanation for the higher 
rotation rates of the galaxies at z ~ 0.5 is that the frac- 
tion of misclassified SO galaxies among visually classified 
elliptical galaxies is larger in intermediate-redshift clus- 
ters than it is locally. If this is indeed the cause of the 
observed effect, then the SO galaxy fraction in the tar- 
get clusters may have been underestimated by 6%. 
Even if this were the strong evolution in the SO 

fraction of clusters continues to be implied between the 
local Universe (where the fraction is ~ 60%) and z w 0.5 
(where for the three clusters studied here the fraction 
is only ~ 22%, even after the listed upward correction). 
A popular explanation for this evolution is that many 
local SO galaxies may have recently transformed from 
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star-forming spiral galaxies. However, this explanation 
is not without problems. Substantial differences would 
be expected between the mean M/L and color of SOs 
and ellipticals, which are not observed (van Dokkum & 
Franx 2001). The contribution of misclassification to the 
observed evolution should therefore remain a topic of in- 
vestigation, in particular since not all misclassified SO 
galaxies would necessarily be recognized from their kine- 
matics alone. 

Mergers can be another reason why there is a decrease 
in the observed average rotation rate between z « 0.5 
and the present. Dissipationless (dry) mergers between 
early-type galaxies can decrease the average rotation rate 
of the population, which is qualitatively consistent with 
the direction of the observed trend. However, it is un- 
clear whether such mergers are numerous enough in the 
high- velocity-dispersion environment of dense clusters to 
reproduce the observed size of the trend. By contrast, 
disk-disk mergers contribute to the comparison of our 
sample to local galaxies only through the so-called "pro- 
genitor bias" . Elliptical galaxies in the local sample that 
formed from such mergers in the last 5 Gyr would not be 
represented in the intermediate-redshift sample, which is 
composed of visually-classified elliptical galaxies. This 
cannot explain the observed trends though, because it 
increases rather than decreases the relative number of 
rapidly rotating galaxies in the local Universe. However, 
such mergers are numerous enough that their progeni- 
tor bias must be taken into account in any quantitative 
attempts to explain the observed trends. 

With large aperture telescopes and multi-slit capabil- 



ities it is now possible to study the resolved kinemat- 
ics of large samples of galaxies at intermediate redshifts. 
As we have shown here, detailed dynamical modeling 
of such data has the potential to provide new insights 
into the evolution of galaxies. This provides important 
advantages over studies that are based solely on global 
quantities. For example, changes in the rotation rates of 
galaxies as a function of redshift can probe their merger 
history. It therefore seems useful to perform more stud- 
ies like the one presented here in the near future, to ex- 
plore wider ranges of redshifts, environments, and galaxy 
types. Such studies would also explicitly constrain the 
contribution to the observed differences in rotation rates 
between local elliptical galaxies and those at intermedi- 
ate redshifts due to potential environmental dependen- 
cies and cosmic variance. 

The models that we have constructed here also yield 
the mass-to-light ratios M/L for all the sample galaxies. 
These are used in Paper II to study the M/L evolution of 
elliptical galaxies as a function of redshift, which in turn 
constrains the formation ages of elliptical galaxies. This 
complements and tests analyses based on Fundamental 
Plane evolution, such as those presented in vDvdM06. 
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